Persistent currents in mesoscopic rings and boundary conformal field theory 
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A tight-binding model of electron dynamics in mesoscopic normal rings is studied using boundary 
conformal field theory. The partition function is calculated in the low energy limit and the persistent 
current generated as a function of an external magnetic flux threading the ring is found. We study the 
cases where there are defects and electron-electron interactions separately. The same temperature 
scaling for the persistent current is found in each case, and the functional form can be fitted, with 
a high degree of accuracy, to experimental data. 
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I. INTRODUCTION 

It is well known that the topology of a condensed mat- 
ter system can strongly influence its excitation spectrum. 
Consider for example the textbook case of a supercon- 
ducting ring. If we write the order parameter in the 
form ip(x) ~ e 1 *^) one can show that there will exist a 
set of topological excitations characterised by the wind- 
ing number 



= — 4 v<i> • <ii 

2?r 



(1) 



where T represents a path encircling the hole in the ring. 
These excitations in general carry a current due to the 
non-zero gradient of the phase field in the ring. This 
effect is very well-studied being fundamentally re- 
sponsible for much of the physics of SQUIDs and related 
devices. 

The prediction of the existence of this "persistent cur- 
rent" is, however, not contingent upon the substrate ma- 
terial being superconducting. Indeed it is apparent that 
the effect should be present in any system with the topol- 
ogy of a ring that permits excitations describable by 
field operators with the C7(l) symmetry alluded to above. 
Persistent currents in normal systems have in fact been 
observed-both in arrays of metal rings O and in indi- 
vidual metal || and semiconducting H rings. From a 
theoretical point of view this is not surprising because 
the persistent current is at heart an Aharonov-Bohm ef- 
fect, independent of the actual nature of the excitations 
in the material. 

This being said, it is quite obvious that the details of 
the supercurrent (such as its amplitude) will be strongly 
influenced by the details of the sample under investiga- 
tion. For example, one expects that the lifetime of the 
excitations carrying the current will be of the order of the 
mean scattering time in the ring, which is heavily depen- 
dent upon the density of scatterers in the material and 
the details of the crystal held. If the associated mean free 
length is much less than the length of the ring L, then 



the effect will never be observed. Rigorous treatments 
of the effects of disorder on the persistent current have 
been performed both numerically || and analytically in 
restricted parameter regimes [^-^|. 

In addition to the effects of impurities, it is suspected 
that the nature of the interaction between the super- 
current excitations can drastically affect the details of 
the resultant current |To|-|T^]. For example, it has been 
shown that at zero temperature and in one dimension, 
the amplitude to backscatter from an impurity renor- 
malizes to zero or to infinity for attractive or repulsive 
inter-quasiparticle interaction respectively |l3| ] . 

In this paper we solve for the partition function of two 
different models; one with arbitrary short-range electron- 
electron interaction strength and no impurities, and one 
with an arbitrary number of impurity sites but no inter- 
action between electrons. These calculations allow us to 
find the finite temperature persistent current developed 
in both cases, which can be directly compared to exper- 
iment H . We find the temperature dependence of both 
models, normalized to their zero temperature values for 
fixed external flux, follows a universal form which agrees 
quantitatively with experiment. 

Our presentation is organised as follows. In section 
II we analyze the case of a defect-free ring threaded 
by a magnetic flux with short-range electron-electron 
Coulomb interactions. We demonstrate that this prob- 
lem can be solved via the use of boundary conformal field 
theoretical techniques. In section III we consider a re- 
lated model, where the Coulomb interactions have been 
removed and an arbitrary number of hopping defects 
added. It is shown that this new model can be mapped to 
one of quasiparticles on a defect-free ring, where the ef- 
fective flux threading the ring is renormalized due to the 
presence of the defects in the bare model. In section IV 
we compare our findings to zero-temperature theoretical 
results and the finite-temperature data of Chandrasckhar 
et.al. H with which we find quantitative agreement. We 
conclude our presentation in section V with a summary 
of our results. 
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II. ANALYSIS OF DEFECT-FREE RING 

A. Introduction of the Non-Interacting Model 

We begin our investigation by considering a tight- 
binding model in one dimension with constant hop- 
ping integrals to and no electron-electron interaction. 
Throughout we shall choose a system of units such that 
h = 1 and fcs = 1. We allow the ring to be threaded 
by a magnetic flux $ = J A^dx = A^L, where A$ is the 
component of the vector potential along the ring. The 
Hamiltonian for this model is 

N 

H = toY J W MN ^]+ii' 3 +h.c.} (2) 

i=i 

with a periodic boundary condition i/jn+i — ipi~ Here 
we define 4>o = 27r<f>/$o with $o = 27rc/e the standard 
quantum of flux and N the number of sites on the ring. 
We assume half- filling (N even); extension to the case of 
N odd is straightforward. 

B. The Continuum Limit and Bosonization 

In this section we rewrite the model in a form that is 
accessible to analysis by conformal field theoretical tech- 
niques. Because conformal symmetry cannot exist in a 
lattice model, the first step is to rewrite (g) as a contin- 
uum model. We begin this procedure by first "opening" 
the ring and introducing an additional site at j = N + 1. 
We then perform a gauge transformation 

_> Jfoi/N^. ( 3 ) 
which allows @) to be rewritten in the form 

N 

Ho=t J2{^+i^+h.c.} (4) 

In this case the periodic boundary conditions are changed 
to 

IpN+l = Klpl (5) 

where k = e^' 1 ^ . We may now rewrite our model in 
terms of left and right moving fermions, namely ip( x ) = 
e lkFX tpL{x) + e~ lkpX ^R(x) 1 where x — ja in the lattice 
model. Taking the continuum limit in the standard way 
gives 

Ho= 2^J dx {^) i-^i> L (x)-{L^R}} (6) 

which is exact up to terms of the order of the lattice 
spacing a. For ease of presentation we choose units such 
that the Fermi velocity Vf = 2toa = 1. It is now possible 
to bosonize @ by introducing a field $ such that 



iPl(x) ~ e -i(fc+*rr») t „ e i(£-a*r») (7) 

on the interval < x < L. Here $ = <&l — is the field 
dual to $ = $l + &r, and r denotes the compactification 
radius of the bosonic field <& ~ $ + 2irr. In the non- 
interacting case under consideration r is taken to be at 
its self-dual point r = ^J= . Note that it is not necessary 
to include Klein factors here as the number of left 
and right moving fermions is conserved in this system. 
Substitution of (^) into the boundary condition (||) gives 

- #(0) = P- , $(L) + $(0)=0 (8) 

We introduce boundaries into the model by folding the 
ring in half. Mathematically this is achieved by introduc- 
ing two bosonic field degrees of freedom defined on the 
half interval 

$ e , o (a0 = ${L - x) ± , Q<x<L/2 (9) 

In terms of the fields $ e ,o(#) the boundary conditions (||) 
become 

* o (0) = ^ , $ e (0)=0 (10) 
and the Hamiltonian becomes 

*-ir*{(50 ,+ (*)'} (ii > 

We have reduced the boundary conditions on the 
fermions (|5|) to simple Dirichlet and Neumann conditions 
on the even and odd bosons while maintaining the free 
nature of the Hamiltonian. 



C. The Role of Conformal Invariance 

Thus far we have translated our original problem of a 
perfect ring threaded by a flux to one of two free bosons 
on the half-interval [0, L/2] with boundary conditions 
(|l0|). The partition function for such a system is imme- 
diately available via the use of boundary conformal field 
theory techniques jlfj . In this section we shall present an 
overview of the relevant conformal machinery and com- 
pute the partition function for our model (|). 

In a conformally invariant theory, a conformal bound- 
ary condition is given by those states which satisfy the 
condition 

( T -T)\ bd .\B)=0 

where T and T denote the holomorphic and anti- 
holomorphic components of the energy momentum tensor 
respectively. This condition corresponds to the physical 
statement that no energy is allowed to flow off of the 
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edges of the system. The general state satisfying the 
above condition is an Ishibashi state fill , 

| M » =£>,iV)®|^iV> 

N 

where the sum is over all the descendents of the heighest 
weight state defined by \fi). In the case of a single free 
boson, with compactification radius r, one can construct 
states which not only satisfy the reflection condition but 
are also eigenvectors of the field or its dual (with eigen- 
value ipo) at the boundary Jl5[ , 

oo ( oo \ 

\D(<p )) = -= e- iWr exp K°' fc )) 

k= — oo I 71=1 J 



M<Po)) =v* E 



a —2irw<po 



exp^ + E44 \ IK )) 



These states correspond to Dirichlet and Neumann states 
respectively. The details of the various operators appear- 
ing here can be found in It is sufficient to note 
that \(w, k)) denotes a Fock vacuum with clq eigenvalue k 
and ao eigenvalue w and a n , a n are the usual annihilation 
operators which make up the bosonic field in canonical 
quantization. 

It is well known that one can compute the partition 
function of a conformal field theory with boundary in 
two equivalent ways 0Q: Z = Tre~ m "- b where H ab 
denotes the Hamiltonian which respects the boundary 
conditions a and b on either end of the system, or through 
a modular transformation, which interchanges space and 
time directions, as Z — (a\e~ LHl3 \b) where Hp denotes 
the Hamiltonian on a circle of radius (3. With this in 
mind, the single free boson has four possible partition 
functions corresponding to the four possible boundary 
conditions DD, NN, DN or ND. These have all been 
computed in jll| and are given by 

Z DD = < D{ Vo )\e-^ H ^D{ip ) > 

17(g) V J 

Z NN = < N(<po)\e-$ H "\N(<p ) > 

l_ q -(A VoM ^ 3 (_ iAvo/ 9 /r;g V2r») 

1 



Z DN = < D(tpo)\e-i H ?\N(<p' ) > 



217(g) 



where Aip = ip — <p' and q = e~ 2 ^P (The length scale 
L has been absorbed in the definition of (3, i.e. to restore 
the length scales (3 — > (3/L). 77(g) and $2,3 (w, q) denote 
the Dedekind eta function and Jacobi theta functions re- 
spectively. Explicitly these are 

00 

v(q) = ? 1/24 II 

n=l 



M™, q) = 2g 1/4 J2 i n{n+1) cos ( 2n + 1 ) t 

n=0 

oo 

i) 3 (!d,g) = l|2^ g" 2 cos 2nto 



n=l 



D. Calculation of the Partition Function 

It is now a straightforward matter to apply these re- 
sults to the free Hamiltonian for even and odd bosons 
with boundary conditions ([To|). Since the two bosons 
are non-interacting the boundary states associated with 
them are tensor products of the single boson boundary 
state. Therefore we can write the relevant partition func- 
tion as 

Z = (U(0)| <g> {N{(f> /2-Kr)\e-^ H e +H ^\N{Q)) ® |D(0)) 
= (D(0)\e-i H e\D(0)) {N{4> Q /2irr)\e-* H e\N{ti)) 



-{<t> /2TT 2 r) 

f] 2 {q) 



-A 



-;q 



(12) 



Notice that we have chosen the boundary state at x = 
L/2 to be <I> = and $ e = 0. One can be convinced that 
this is correct by starting with the condition (||) where 
the sites N + 1 and 1 are replaced by N/2 + 1 and N/2 
respectively and taking k = 1. This effectively sets the 
flux in (|l^) to zero. 



E. Inclusion of Coulomb Interactions 

This completes the analysis of our original toy model 
(||). Let us now consider what happens when we include 
short-range Coulomb interactions. We repeat the proce- 
dure followed in the previous sections with the Hamilto- 
nian 

N 



+ tfoXywcfl^o' + i)<Ki + 1) (13) 



It is well known that the only effect of the Coulomb 
term is to renormalize both the radius of compactification 



V47T V t 7T 

and the Fermi velocity 

Up 
2vrf 



Vf = 2t$a 1 + 



(14) 



(15) 
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In order to accomodate this we choose new units such 
that Vf = 1. It is then straightforward to use (|l^ ) to cal- 
culate the free energy of the system and then the persis- 
tent current as a function of the compactification radius 
r{Un). Explicitly this is 

edlnZ , s 

p O(p 

We show in fig.l the results of this calculation. Varying 
the parameter r does not qualitatively affect the shape 
oiJ{<h,P). 



J 




FIG. 1. Persistent current in a defect-free ring for arbi- 
trary Coulomb interaction strength (units for (5 are chosen to 
highlight the transition from the high temperature regime to 
the low temperature limit (corresponding here to /3 ~ 1) and 
units for J are arbitrary). Note that J is periodic in <j>o with 
period 2-k. 



III. INCLUSION OF HOPPING DEFECTS 

We now demonstrate how this analysis may be ex- 
tended to deal with the presence of hopping defects. The 
Hamiltonian for a ring with m defects reads, after an 
appropriate gauge transformation, 



H 



(17) 



Here the set of sites {^4} labels the various defect links 
denotes the strengths of the defects. 



and k, 



Pie 



i<j>o/r. 



It is a straightforward matter to diagonalize (17); 



Ho 



fc:/(fc)=0 



(18) 



The dispersion relation here is = 2iocos(fc), and the 
restriction to f(k) = implements the quantization con- 
dition. Note that in all that follows we will only be in- 
terested in energies very close to the Fermi surface. We 



shall assume that if there are any localized states, their 
energies are far from Ef and their contribution to the 
low-energy physics can be neglected. 

As an explicit example of the quantization condition, 
one finds for the case of a single defect f(k) — reduces 
to M: 



f 1 (k) = p 2 sm(k(N-l)) 
+ 2pcos((/>o) sin(/c) 



i(k(N + 1)) = (19) 



In the low-energy limit under investigation, for an arbi- 
trary number of defects, the quantization condition may 
be linearized near the Fermi points k = zt^. Performing 
this procedure, one finds that the effects of the defects 
and flux may be subsumed into a shift in the fermionic 
spectrum given by 



k = ±- 



2tt a{{K},(j> ) 
N U+ N 



n€7L. 



(20) 



For the single defect case ( |T9| ) one can show 20 that a 
is given explicitly by 



2 P 



1 



■ COS <po 



(21) 



In general for m defects a is found by solving a secular 
determinant of order 2™\ 

Introducing new quasi-particle fermionic operators 



e N 



-TT'iJ 



(22) 



into our Hamiltonian (fi~8[) we find that 



N 

H = -hJ2{ eia/N ^ti + 1)^0') + h - c ) (23) 

which is the original toy model (||) with an effective flux 
given by a instead of <po- This demonstrates that the 
low energy limit of the multi-defect case can be reduced 
to one of free quasi-particles on a ring threaded by an 
effective flux a({n}, O ). At this point one can apply the 
analysis of the previous section to obtain the partition 
function and therefore the persistent current. Shown in 
fig. 2 is the persistent current in the simplest case where 
we have only one defect. We find that the flux depen- 
dence of J changes dramatically as p departs from 1, 
evolving from the sawtooth shape seen in the no-defect 
case to a sinusoidal variation. 



IV. COMPARISON OF TEMPERATURE 
DEPENDENCE OF PERSISTENT CURRENT TO 
EXPERIMENT 

We have calculated the partition functions for both 
models under consideration. This allows the temperature 
dependence of the persistent current to be extracted for 
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FIG. 2. Persistent current developed in a ring with no 
electron-electron interaction and one hopping impurity with 
strength p = 6. Units for J and (3 are the same as those 
defined in figure 1. 



model on the region [— 7r..7r]. It is interesting to note 
that the form of ( |26| ) is not a strong function of cv, its 
shape nearly entirely dependent on the value of T\ (see 
fig]^)- This means that regardless of the model chosen, 
this one parameter (which can be varied arbitrarily in 
both cases) determines the temperature dependence of 
the normalized persistent current for fixed </>o, regardless 
of the resultant value of a. 
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each. We may perform both calculations in parallel by 
retaining all factors of r and a in our expressions, keeping 
in mind that these cannot be varied independently. For 
the defect-free model, a = 4>o and r is arbitrary. For the 
non-interacting model, a is given by solving for the roots 
of the quantization condition f(k) = (Eq) linearized 
near the Fermi points and r = 1/\/4tt. 

With this caveat, the free energy can be written in 
general as 



F = - 1 -mz 



2tt 3 t 



i(lntf 3 (0;e^ 2 ) 



ia/3 



l n tf 3 (-^ ;e -^A )_ 2 l n?? (<z) 



with subsequent persistent current 
J 



dF _ da OF 
OA d<f) da 



whose zero temperature limit is simply 

ea da 



(24) 



(25) 



If we normalize the persistent current to this value we 
find that 



a 3 \2n 2 T ,e 



(26) 



where we have introduced the energy scale 



Tr 1 



Note that though a can be a complicated function of 
</>o and the defects, it is bounded, being defined in our 



FIG. 3. AJ max is the maximum variation in J/ Jo for 
a G [— 7r..7r], here given as a function of Ti/T. We see that 
AJmax « 1 for all Ti/T, demonstrating the weak depen- 
dence of J I Jo on a. 



We now compare the form derived for J/ Jo to the ex- 
perimental results of Chandrasekhar et.al. j^] obtained 
from mesoscopic gold rings. In these experiments the 
authors present data for the maximal normalized persis- 
tent current, defined to be the solution of (26) with (j) 
chosen so as to maximize the zero to peak amplitude of 
the signal. Shown in fig. (Q) are the results of a fit to 
the data. The circles represent values gathered from a 
1.4/im x 2.6/iTO loop and the squares the response from 
a 2.4/7777 ring. We find that the value T x = 0.09K fits 
both cases extremely well. This finding lends credence 
to our claim that the details of the defects do not af- 
fect the temperature dependence of J/ J . This is simply 
because the two samples considered in the experiments 
must have had differing defect structures (corresponding 
in our language to different shifts a) yet the experiments 
show the same temperature dependence. 



V. CONCLUSIONS 

Using boundary conformal field theoretical techniques 
we have computed the partition function of a tight- 
binding model of electron dynamics in a mesoscopic nor- 
mal ring. Defects and electron-electron interactions were 
introduced separately. In the case where defects are 
introduced, we mapped the problem onto a defect-free 
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FIG. 4. Fit to the data of Chandrasekhar et.al. 
Tl = 0.09K. 



with 



model with an effective flux which depends on the details 
of the original defects. The second case, where electron- 
electron interactions are included within the defect-free 
model, was solved by noting that such interactions serve 
only to alter the compactiheation radius of the bosons 
and the effective Fermi velocity. 

Using the partition functions we computed the per- 
sistent current for both cases as a function of flux and 
temperature. The model with defects generated a per- 
sistent current that was found to change its flux depen- 
dence from a sawtooth to a sinusoidal shape even for 
one weakly scattering impurity. Futhcrmore, the am- 
plitude of the oscillations decreased as the strength of 
the defect increased. However, the inclusion of electron- 
electron interactions in the defect-free case produces a 
different behaviour: the sawtooth nature is not altered 
and the amplitude of the current decreases only very 
weakly as the interaction strength is increased. Although 
the functional dependence on flux differs in the two cases 
described here, their normalized currents have identical 
temperature dependence. There is one single free parame- 
ter in the expression for the current which sets the scale in 
which temperature is measured. This suggests that such 
a functional dependence on temperature is generic in the 
sense that the microscopic details of the interactions and 
impurities do not affect it when taken separately. 

We compared the temperature dependence with ex- 
periment |3|] and with only one fitting parameter (the 
energy scale) found excellent agreement. In such an ex- 
perimental setup one expects that both electron-electron 
interactions and defects play a significant role, yet the 
predictions of our simple model were able to capture, 
quantitatively, those results. This seems to indicate that 
the results obtained for the normalized persistent cur- 
rent are not sensitive to any interplay between defects 
and electron-electron interactions in the bare model-that 
is, including defects and electron-electron interactions si- 
multaneously should lead to analogous predictions. 
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